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ABSTRACT 

We combine data from SDSS and the GALEX and Spitzer observatories to create a sample of galaxies ob- 
served homogeneously from the UV to the Far-IR. This sample, consisting of ^460 galaxies observed spectro- 
scopically by SDSS provides a multiwavelength (0. 15-24 /jm) view of obscured and unobscured star formation 
in nearby (z < 0.3) galaxies with SFRs ranging from 0.01 to 100 M Q yr" 1 . We calculate a robust dust measure 
from the infrared to UV ratio (IRX) and explore the influence of star formation history (SFH) on the dust-UV 
color relation (i.e. the IRX-/? relation). We find that the UV colors of galaxies are only weakly dependent 
on their SFH as measured by the 4000A break. However, we find that the contributions of dust and SFH are 
distinguishable when colors at widely separated wavelengths (e.g. 0.23-3. 6pm) are introduced. We show this 
explicitly by recasting the IRX-/3 relation as a more general IRX-SFH-color relation, which we examine in 
different projections. We also determine simple fits to this relation. 

Subject headings: galaxies:evolution — dust:extinction — galaxies :ISM — ultraviolet:galaxies — in- 
frared:galaxies 



1. INTRODUCTION 

The color of a galaxy is determined primarily by its star 
formation history (SFH) and the amount of dust attenuation 
present, with significant additional contributions from metal- 
licity and dust geometry. An empirical relation between SFH, 
dust attenuation, and color is therefore a useful constraint on 
models of galaxy formation, and can aid in the interpretation 
of high redshift galaxy observations where measurements are 
difficult. 

Here we investigate such an empirical relation, using a rep- 
resentative sample of galaxies observed from the ultraviolet to 
the infrared by GALEX, SDSS, and Spitzer. The long wave- 
length coverage allows us to construct simple but robust mea- 
sures of dust attenuation that are relatively free of a depen- 
dence on SFH. Similarly, the SDSS spectroscopy allows us to 
measure SFH diagnostics - we use a 4000A break measure - 
that are largely free of a dependence on dust attenuation. 

Earlier studies of the effect of dust on galaxy colors have 
focused on the ultraviolet (UV, AA ~ 1400-2500) colors of 
galaxies. A primary reason for this is that the intrinsic, un- 
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derlying UV color (before attenuation) is relatively insen- 
sitive to the SFH of the galaxy, when compared to the ef- 
fects of dust attenuation. Also, studies of the highest redshift 
galaxies are often restricted to the restframe UV so that es- 
timate s of dust atten uation must be made using UV colors. 
Calze ttiet all (11994) derived the effective attenuation proper- 
ties of dust in star-bursting galaxies by comparing the total 
dust absorption - measured from the ratio of dust emission 
to UV emission - to the change in UV color of a sample of 
starforming galaxies (the so-called IRX-/3 relation). This lo- 
cally derived relation has been used extensively to determine 
the dust attenuati on in galaxies at hig her redshift using their 
UV colors ( e.g. iMe urer et a ll dl999l)J. However, more re- 
cent studies ( Bel l et alJl2002t LSeibert et alJl2005t iBuatetalJ 
2005) have shown that there is significant scatter in this re- 
lation, especiall y when less r apidly starforming galaxies are 
included. Kon g et all y004) show, using stellar population 
synthesis models with significant recent bursts, that both SFH 
and dust are expected to affect the UV color of normal galax- 
ies since at a constant value of dust attenuation galaxies with 
an older stellar population should appear redder in the UV be- 
cause of their redder intrinsic UV spectra. They thus propose 
that a measure of SFH - specifically the 4000A break - can 
be used to explain the scatter in the LRX-/3 relation. 

2. DATA 

We have constructed a sample of galaxies observed from 
the UV to the mid/far-infrared. Our primary data set is 
the sample of local (z < 0.3) galaxi es observed spect r oscop- 
icall y by SDSS and a nalyze d by iKauffmann et all (2003) 
and iBrinchmann et all (120041) (hereafter SDSS/MPA galax- 
ies). The UV data is taken from pipeline processed GALEX 
observations of the Lockman Hole, with exposure times of ~ 
1.5 ks. The pipeline-produced GALEX catalogs are searched 
for objects within 3" of the SDSS/MPA galaxy locations, and 
the nearest object is taken as the match to the SDSS galaxy. 
UV flux measurements are made in elliptical Rron apertures. 
In the optical we use SDSS petrosian magnitudes. The in- 
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FIG. 1. — IRX vs. color of the sample for different bins of D„(4000). Left: f — n color. Lines show the relations of Meurer et al. 1 1999) (dashed) and 
Seibert et al. 1 2005) (solid) where the conversions from Af m and (3 are from Seibert et al. 1 2005). Middle: The long-baseline n — 3.6/im color. Dashed line gives 
the fit from the 3rd row of TablelTlfor D„(4000) = 1 .25. Right: The distribution of D„(4()00)for the sample, with dashed lines and color bar showing the quintiles 
used to color code the sample galaxies. Error bars show the median measurement errors (excluding systematic effects) for the entire sample (black) and the 
lowest and highest quintiles of D„(4000) (purple and red respectively). Note the clear separation of galaxies with different D„(4000) in the middle panel, while 
for the / - n color (left) such a separation is marginal. 



frared data is provided by the SWIRE Spitz er observations 
of the Lockman Hole (jjLonsdale et al. 2003). We have per- 
formed aperture photometry in the SWIRE team processed 1 
3.6 through 7.8/im IRAC images and 24/im MIPS images 
at the location of each of the SDSS/MPA galaxies, using a 7" 
radius aperture (12" at 24 /im). The fluxes are then aperture 
corrected to total magnitudes. Systematic errors in IR flux 
due to calibration uncertainty, aperture corrections, and the 
resolved nature of many of the sources amount to ~ 30%. 

The resulting UV through 3.6/im magnitudes are K- 
co rrected to z = 0.1 (e. g. u, 0A g, etc.) using the method 
of iBlanton et al.l J2003I) . At longer wavelengths dust emis- 
sion becomes more important than stellar emission, and we 
use a different method to '^-corr e ct' the data: we choose the 
best fitting redshifted iDale et alJ (120011) model IR SED, on 
the basis of the observed 8 to 24/im flux ratio. This SED is 
then normalized using the measured 24/im flux, and the in- 
tegrated far-infrared (8-1000/i m) dust luminositie s (Lj„ sf ) are 
derived. Note that the different Dale et all (120011) SEDs have 
Li4fim/Ldust ratios that are different by a factor of up to five. 
We have checked that our results woul d not change signifi- 
cant ly if we use the model SEDs of Devriendt et al. ( 1999) 
(see Papovich & Bell (2002) for a detailed discussion of pre- 
dicting IR luminosities from Spitzer data). 

Our final selection of galaxies consists of those SDSS/MPA 
galaxies with detections in the Far-UV (/, A = 1516A) and 
Near-UV (n, A = 2267 A) through 24/im bands. This is a 
total of 467 of the 645 SDSS/MPA galaxies within the ~ 9 
deg 2 of SDSS/GALEX/Spitzer overlap. These galaxies have 
a stellar mass range of ~ 10 7 - 10 1L5 M and a SFR range 
of ~ .01 - 100 M^yr" 1 , as d e termin ed bv [kauffman n et aTl 
( 2003 ) and Brinch mann et al.l (120041) from the optical spec- 
tra and photometry. To simplify the present analysis we do 
not consider upper limits or selection effects except to note 
that we are biased against galaxies with very low Lf uv or 
Ldmt- Nevertheless, we recover a significant number of galax- 
ies which appear 'old' and elliptical but have low levels of UV 
and IR emission. 

1 http://swire.ipac.caltech.edu/swire/astronomers.html for a description of 
the SWIRE image processing 



3. ANALYSIS 

3.1. Dust Indicators 

One of the primary motivations for compiling the sample 
of galaxies described above was to construct a robust and 
model independent measu re of dust attenuation , the so-called 
infrared excess (IRX, see iGordon et alJ (2000) for a discus- 
sion of the relation of infrared excess to UV attenuation). We 
adopt the definition IRX = log(Lj MJ , /L f uv ) where Ldust is the 8- 
1000/im dust luminosity as determined above and Lj m , = vh v 
is the luminosity in the 01 f band (v w c/1400A). A sec- 
ond availa ble dust indicator is the Ha to H/3 decrement (see 
iKennicutH (119981) and references therein), measured from the 
SDSS spectra. This dust measure is only well defined for 
those galaxies with strong emission lines, and cannot be eas- 
ily compared to global galaxy measures due to spectroscopic 
aperture effects. We do not consider it in the present study. 

3.2. SFH Indicators 

We take the str ength of the 4000A break in galaxies as an 
indicator of SFH dBaloeh et al.H999tlBrinchmann et al.l2004l 
Mac Arthur 2005|. This has been measured from the dered- 
shifted SDSS spectra by iKauffmann et all (120031) using the 
ratio of the flux in two narrow bands (AA = 100) centered at 
4050A and 3900A. This narrowband color, D„(4000), is less 
sensitive to reddening by dust than broadband col ors. It is not, 
howe ver, completely insensitive to dust effects (MacArthur 
120051) . Note that D„(4000) is only measured within the 3" 
SDSS aperture, which can cause an overestimate of the in- 
tegrated D„(4000) for galaxies with moderate bulge to disk 
ratios. We do not consider in detail the relation between 
D„(4000) and more physical measures of galaxy age or SFH 
(e.g. the ratio of current to past averaged star formation rate 
or the specific star formation rate), which is metallicity de- 
pendent and requires population synthesis modeling. 

4. RESULTS: THE DUST-SFH-COLOR RELATION 

4.1. IRX vs. Color Binned by D„(4000) 

Models of galaxy spectra have suggested that the UV color 
of galaxies can be decomposed into contributions from dust 
and SFH, thou gh the eff ect of dust is the dominant contri- 
bution dKong et alJl200l . In Figure [2 we show the relation 
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between IRX and color for different ranges of D„(4000). We 
find a marginal dependence of the IRX vs. f—n relation (anal- 
ogous to the YKX-f3 relation) on D„(4000). Only the very old- 
est galaxies (red points) have systematically redder f — n for 
the same value of IRX, but also have larger errors. For galax- 
ies with D„(4000)< 1.8 there appears to be little dependence 
of the scatter on ga laxy SFH. This is similar to the result of 
ISeibert et alJ (12005), who use the n-K color as a measure of 
SFH (but see below for the effect of dust on the «-NIR color). 

This result can be reconciled with the models of Kon g et aTl 
(2004) in several ways. First, some stars contributing to the 
dust heating (and therefore the IR emission) may not be con- 
tributing to the UV emission or, conversely, some stars con- 
tributing to the UV emission are not contributing to the IR 
emission (i.e. a decoupling of IR and UV emission). This 
may occur for var i ous reasons including star -dust geometry 
(e.g. lCalzetti et al.l(l2005l) : lThiiker et al J l l2005l) ') and contribu- 
tions to dust heating from older stellar populations. Related to 
this, D„(4000) may not be measuring the relevant timescale 
for changes in the UV spectrum. The SFH diagnostic that 
is likely to be most relevant to changes in the UV spectral 
shape is the ratio of current SFR to the SFR averaged over the 
last 100 Myr dCalzetti et al.ll2005l) . which is not probed well 
by D„(4000). Second, weak AGN may affect the UV colors. 
Third, the entire range of f — n color is < 2 mag, compared 
to a median error of ~ 0.25 mag. The scatter due to mea- 
surement erro r may obscure a tre nd in the scatter with SFH. 
The results of Kon g et aTl (120041) were based on a library of 
95000 model spectra, which may not be well represented by 
our sample. These and other possibilities will be investigated 
in future work. 

For n — 3.6/im color it is easy to s ee the effect of SFH that 
was predicted bv lKong et alJ |2004) for the f—n color. At a 
given low D„(4000) (i.e. younger mean stellar age, purple 
points) the relation between dust and color is clear, and has 
low scatter. This relation is closely related t o the eff ective 
attenuation curve of the d ust in these galaxies ( Calze tti et alJ 
119941 iMeurer et alJ ll999). For intermediate D„ (4000) (i.e. 
intermediate mean stellar age, orange or green points) the re- 
lation between dust and color persists with low scatter, but the 
entire relation is shifted to redder color. This is presumably 
because of the redder intrinsic spectrum of an older stellar 
population, on top of which the effect of dust attenuation on 
color persists relatively unchanged. For n— 3.6/xm the ratio of 
the median error to the range in color is significantly smaller 
than for f—n, making trends with color easier to see. 

At the largest D„(4000) the scatter in the IRX-color re- 
lation increases - this may be due to several causes. First, 
there are larger errors in IRX (and color) for these galax- 
ies, which have relatively little star formation and are sys- 
tematically dimmer in / and at 24/xm. Second, the effects 
of metallicity on D„(4000) become more pronounced at high 
D„(4000). Third, some UV emission may be due to evolved 
populations (e.g. BHB stars). Finally, these galaxies may well 
host AGN that affect the IR and/or UV emission, changing 
IRX. 

4.2. D„(4000) vs. Color Binned by IRX 

In Figure|2]we present another projection of the dust-SFH- 
color relation. Here we show D„(4000) as a function of 
«-3.6/imcolor for different ranges of IRX. IRX appears to be 
well determined using just «-3.6/im and D„(4000). The rela- 
tion between D„(4000) and n — 3.6/xm shifts to redder color 
for galaxies with more dust, while the slope of the relation 




FIG. 2. — D„(4000) vs. n — 3.6/^m color for the sample, with color coding 
by IRX. The colored lines show a fit to the relation (3rd row from Table fTT. 
where each line represents the fitted relation for the median value of IRX 
in the bin of the corresponding color. Labels indicated the type of galaxies 
typically occupying different regions of the diagram (SF-star formation; RS- 
red sequence). Error bars show the median measurement errors (excluding 
systematic effects) for the entire sample (black) and the lowest and highest 
bins of IRX (purple and red respectively). Inset, Lower Right: Distribution 
of IRX for the sample. Dashed lines and color bars show the bins used to 
color code the the sample. 

remains nearly constant. Galaxies in the lower right of this 
plot are predominantly dusty star-forming galaxies. Galaxies 
in the lower left are typically blue star-forming galaxies or 
dwarfs/irregulars. Galaxies in the upper right of Figure|5]are 
red-sequence galaxies composed primarily of old stars, with 
very low levels of UV and IR emission. 

4.3. Fits to the Relation 

We have conducted quantitative parametric fits to the rela- 
tion between dust attenuation, SFH, and the «-3.6/im color. 
Considering a simple effective attenuation law we are moti- 
vated to consider polynomial fits to the color (c), D„(4000), 
and an expression of IRX that is linearly related to color 
A IRX = 2.51og(10 /sz - logBC * + 1) where 5C* = 1.68 is the bolo- 
metric correction to the UV and the bolo metric correctio n 
to the IR ha s been made in ^3 ( see |Meurer et alJ ( 119991) : 
iGordon et all ( l2000l) : ISeibert et all ( 120051) for a discussion of 
the relation between IRX and the true UV attenuation Apuy). 
We treat the color and D„(4000) as independent variables 
(since the errors are much smaller than for IRX). We assume 
either a linear or quadratic form for the relation between color 
and D„(4000). The results are given in Tabled Examina- 
tion of the residuals shows that a cross term ((D„(4000)) * c) 
is required, and fits including such a term are also given in 
Table □ The fit from the 3rd row of Table fTI is overplotted in 
Figure [2] and the residuals in IRX as a function of color and 
D„(4000) are shown in Figure [3] Figure [5] shows that the fit 
is poorer for redder colors and larger D„(4000). This may be 
due to the effects listed in 34.1 I f or the » -3.6/ im color, as well 
as the effect of aperture on D„(4000) ( 33.2b . At bluer colors 
(n-3.6/im< 4) and smaller D„(4000) the residuals are lower, 
~ 0.4 rms in Ajrx- Table^also includes fits using the n-r 
color since this is more easily measured for a large sample o f 
GALEX observed SDSS galaxies (e.g. ISalim et alJ f2005)1. 
The behavior is similar to the 7i-3.6/im color though the 
residuals are larger. 
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TABLE 1 

Fits to the Dust-SFH-color Relation of the form 11 
Airx = A + Bjt+O 2 +D;y + Exy WHERE x = (D„(4000)-1.25) 
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FIG. 3. — Residuals in IRX from a simple functional fit of the dust-SFH- 
color relation (3rd row of Table HI. Dashed lines show the smoothed rms 
dispersion of the residuals. Left: Plotted against the n-3.6^tm color. The 
color of each point indicates its D„(40()0) quintile. Right: Plotted against 
D„(4000). The color of each point indicates its n — 3.6/xm color quintile. 

We have used a simple treatment of attenuation to motivate 
our fits that does not self consistently consider, e.g., the heat- 
ing of dust by old stars. Also, D„(4000) is subject to aperture 
effects that may affect the fits. It is thus difficult to interpret 
the fit coefficients as physical parameters. The fits in Table^] 
are referenced to D„(4000)= 1 .25, c = 2.0 because this region 



is typical of 'blue sequence' galaxies. The constant A term 
thus gives the typical Airx for such galaxies. The B term (with 
contributions from the C and E terms when present) gives the 
trend in Airx with D„(4000), and is negative since galaxies 
that are intrinsically redder have less attenuation for a given 
color. Similarly, D gives the trend in Airx with color, and 
is positive as a redder galaxy at a given D„(4000) has more 
attenuation. 

5. CONCLUSIONS 

While we recover only a weak SFH dependence for the 
short wavelength-baseline UV color, we have shown that long 
wavelength-baseline colors - specifically n-3.6/im - can 
be decomposed into contributions from dust and SFH with 
low scatter. This is the dust-SFH-color relation. Such a de- 
composition is possible due to the use of a relatively dust- 
insensitive (though perhaps metallicity dependent) SFH indi- 
cator in combination with a robust measure of the dust atten- 
uation in galaxies. At large D„(4000) the relation between 
dust, SFH, and color is more scattered, suggesting that an ad- 
ditional parameter may be necessary to explain the color of 
such galaxies, or that our SFH and/or dust indicators become 
less reliable here. Deep Spitzer IRAC data, when combined 
with optical data, is ideally suited to measuring the n-3.6/im 
color for galaxies at z = - 2, allowi ng similar analyse s for 
galaxies at much e arlier epochs (e.g. iKriek et afl J2006I) and 
iReddv et alJESol ). 

The MPA/JHU collaboration for SDSS studies has very 
generously made their catalogs publicly available. The pub- 
licly available Spitzer data obtained and reduced by the 
SWIRE team have been essential to this work. We gratefully 
acknowledge NASA's support for construction, operation, and 
science analysis for the GALEX mission, developed in coop- 
eration with the Centre National d' Etudes Spatiale of France 
and the Korean Ministry of Science and Technology. BDJ was 
supported by NASA GSRP Grant NNG05GO43H. 
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